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We establish the connection between the presence of a glass phase and the appearance of a 
Coulomb gap in disordered materials with strongly interacting electrons. Treating multiparticle 
correlations in a systematic way, we show that in the case of strong disorder a continuous glass 
transition takes place whose Landau expansion is identical to that of the Sherrington-Kirkpatrick 
spin glass. We show that the marginal stability of the glass phase controls the physics of these 
(-H ' systems: it results in slow dynamics and leads to the formation of a Coulomb gap. 
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The relation between the slow dynamics of Coulomb glasses, disordered materials with strong electron-electron 

, repulsions, and the appearance of a soft "Coulomb" gap in their density of states (DOS) has been a mystery for a 

long time. The strong effect of interactions on the DOS was first noticed by Pollak yj and Srinivasan |2j]. Efros and 
Shklovskii J3J predicted the Coulomb gap as resulting from the long-range Coulomb interactions between localized 
electrons in semiconductors, and leading to a crossover in the temperature dependence of the conductivity from 

1^ Mott's law ln(p) - (Tm/T) 1 / 4 to ln(p) - (Tes/T) 1 / 2 at low temperatures 0. The Efros- Shklovskii law for the 
conductivity was verified in semiconductors, alloys and granular metals, and recently the gap itself was directly 
observed in semiconductors [a, li| ■ The presence of disorder frustrates the Coulomb interactions and leads to glassy 
behavior in such materials, as predicted theoretically long ago pj. The first evidence for glassy phenomena came from 
the slow relaxation of charge injected into compensated semiconductors Jjjj. Later, Ovadyahu's group established 
that the slow dynamics in Indium-oxides is indeed due to electron-electron interactions [9j- Very recently the same 
group has demonstrated memory and aging effects similar to those observed in spin glasses | lOj . 

Despite the experimental progress a thorough understanding of the glass phase is still missing. The source of the 
difficulty is that in order to describe glassy phenomena one needs to take electron correlations into account, while 
the approach by Efros and Shklovskii is essentially a single particle theory. The necessity to include correlations has 
£> . become clear from several recent numerical studies of the off-equilibrium dynamics of Coulomb glasses [ll|. Further, 
an increasing number of experiments suggests that a quantitative description of the hopping conductivity should take 
multiparticle processes and correlations into account |12J . One also needs to go beyond a single particle theory in 

^1J ' order to understand the relation between the Coulomb gap and the glass transition. The mean field solution of a 
model of uniformly interacting electrons in a disordered medium indicates that the glass transition and the formation 

r^- of a pseudogap in the DOS are driven by the same mechanism, and a similar relation has been conjectured for the 

<0 Coulomb glass J13J. 

The goal of this paper is to develop a formalism that accounts for the correlations between the electrons in a realistic 
model for Coulomb glasses in 3D. Our approach is based on the locator approximation developed for spin glasses in 
Refs. |l4J ll<4 llfij. The idea is to map the original lattice problem into an effective single-site problem that encodes 
correlations by the distribution of a fluctuating local field, which gives exact results for infinite range models. In the 
limit of strong disorder, the Coulomb interactions are essentially unscreened, so that the effective number of neighbors 
is large and the locator approximation is parametrically well justified. In this regime, we find a replica symmetry 
breaking glass transition, which belongs to the same universality class as the transition in the Sherrington-Kirkpatrick 
^ . (SK) spin glass |17| . Below the transition, the phase space divides into an exponential number of metastable states 
and ergodicity is broken. Like any generic glass, this electronic system freezes into one of many marginally stable 
states since the latter are the most abundant (the number of states increasing exponentially with decreasing stability). 
Marginally stable states are characterized by the presence of soft modes that lead to the slow relaxation dynamics 
observed in experiments. Above T c , the DOS does not display any particular signature of the Coulomb interactions. 
We show that the Coulomb gap forms only below T c where it emerges as a direct consequence of marginal stability. 
Finally, we derive an asymptotic expression for the DOS at very low temperatures. 

We consider the classical model |4| for strongly localized electrons occupying a fraction < K < 1 of a given set of 
impurity sites i, 

H=-^2n i Jijn j + ^2n i (e i -n K ), (1) 

i^j i 

where n, € {0, 1} is the occupation number of the site i. For simplicity, we take them to be arranged on a cubic lattice 
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FIG. 1: For long-range interactions the self-energy E can be approximated by a local operator. The full propagator is obtained 
as a simple geometric series. 

with lattice spacing I = 1. The unscreened Coulomb interactions are described by Jij — 1/rij in units where e 2 /£ = 1. 
The €i's denote random on-site energies, and fj, K is the chemical potential. We restrict ourselves to the case K = 1/2 
where the particle-hole symmetry implies [i 1 / 2 = 0, and suggests to introduce spin variables Si = rii — 1/2 = ±1/2. 
Further, we assume a Gaussian distribution of width W for the on-site energies e^. Their randomness emulates the 
effect of the disorder in the site positions which is present in all physical electron glasses and generates rather large 
site-to-site fluctuations of the Coulomb potential. We focus on the limit of strong disorder, W 3> 1 and dimension 
D = 3. In this case screening is suppressed on short scales and the interactions remain long-range, which justifies the 
use of the locator approximation. We will see that at low temperatures the self-generated disorder outweighs W. We 
thus expect our results to be universal in that regime. Furthermore, this observation makes us believe that at low 
temperatures the locator approximation can be justified even in the case of weak disorder. 

In the case of long-range interactions, diagrammatic expansions can be efficiently resummed since the large number 
of effective neighbors allows one to approximate the self-energy by an average local term, see Fig.^ This observation 
further suggests to replace the interactions of a given spin with its environment by an effective local field described by 
couplings of the spin to its own replicas. This reduces the model to a single-site probl em, translating the complexity 
of the environment into a non-trivial replica structure of the one-site Hamiltonian [13J, I1(t| , 



f3H ({s a }) = -^s a B ab s b +/3y% a e . 



(2) 



a , b 



We note that this procedure resembles the way in which the SK-model is transformed into a one-site problem. 
Indeed, the locator approximation applied to the SK-model sums all tree-like diagrams with doubled interaction lines 
which becomes exact in the large N limit. 

To make contact between the Hamiltonians |JQ and (J2J, we require that they both yield the same single-site spin 
correlation functions, 
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We have averaged over the random energies e^, and as usual, the number of replicas n is implicitly assumed to tend 
to zero in the end of calculations. X denotes a n x n block matrix with all entries equal to 1, and we have defined 
B — B — f3 2 W 2 l. In Eqs. © we approximated the full propagator for either model as a simple geometric series with 
a local self-energy insertion £ f,<$ij, as motivated above. Since the mapping is to preserve correlations, the self-energy 
has to be the same for both models [16(. From IJ2J), we obtain the free energy 
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(4) 



where U(B) has to be determined such that the saddle point equations with respect to B yield back Eqs. J3J. Up to 
a function of temperature, we find 



U{B) = tr 



\n(B - E) - -^Tr In ((3 J - (3 2 W 2 1 - E) 



(5) 



where tr denotes the trace in replica space. We emphasize that in this expression the self-energy E has to be 
considered as an implicit function of B as defined through Eq. @ . In the following we shall need spatial traces like 
g n {x) = V~ l Tr , j\ x)n which we evaluate in Fourier space g n (x) = J tjMts , 0J *>„ with Jk = 4n/k 2 at small k. We 
assume some cut-off procedure that regularizes the small scale physics so that J . = 1 and J, Jk = 0. For x 3> /3 we 
obtain <7„(x) « x _ ™[l + C„(/3/x) 3 / 2 ] where C\ = 2y/n and C2 = by/ir. 

Let us first discuss the replica symmetric (RS) solution of Eqs. (01 for which we assume £ a & = — ^oS a b + EjZ, 
and B a b = —B 5 a b + Bjl. For If > 1, we find E ~ \/2ttPW, suggesting the interpretation of Eq -1 as the 
fraction of thermally active sites (for T c < T <C W) . The distribution of local fields obtained from this RS solution 
matches remarkably well numerical simulation data for finite temperatures and W < 1, even though the locator 
approximation is difficult to justify in this parameter regime |l8j . A depletion of sites in small fields is found due to 
strong correlations in this "Coulomb plasma" . However, a closer analysis reveals that there is no true pseudogap on 
the replica symmetric level, the depletion disappearing completely for strong disorder. This is also reflected in the 
charge susceptibility Xrs = 0(1/4 — (s a Sb) a m) ~ 0/^o which tends to a finite constant (~ 1/W) within this solution. 
The genuine Coulomb gap is formed only when the replica symmetry is broken. For W ^$> 1, the RS solution indeed 
exhibits an instability when the condition 

e -y 2 /2(W 2 +B x /p 2 ) 1 

■ d V 2 ^ 2 +^// 32 ) [2cosM ^ /2)J4 (6) 

= [^(EoJ-fc-^Eo)]" 1 « M 3 So)- 1/2 , 

is met, from which we extract the critical temperature T c w W^ 1 ^ 2 /[6(2/tt) 1 / 4 ] <C 1 ^C W. We emphasize that the 
difference 5 1 " 2 (Eo) — ff^^^C^o) is controlled by the contribution from large scales, 1/fc ~ VW, which justifies our 
assumption of a large number of effective neighbors. 

The instability © signals a continuous glass transition with full replica symmetry breaking. We may analyze it 
further by expanding the free energy with respect to the replicon mode SB (with SB aa = and 5B2 — 0) around the 
RS solution, 



nf38F = n 



w z i 2 



tr(~rSB 2 + c 2 5B 3 ) + c 3 ^ SB, 



(7) 



where r = 1 — T/T c . This shows that the glass transition in Coulomb glasses belongs to the same universality class 
as the one in the SK-model. Hence, many results about the critical behavior known for infinite range spin glasses |l9j 
should be directly applicable to the present case. This might be interesting in particular for the aging and memory 
effects observed in experiments 10], even though the dynamics of spin glasses obey slightly different rules than in 
Coulomb glasses. Vice versa, electron glasses present an appealing testing ground for many theoretical ideas developed 
in the context of the SK-model. 

We now turn to a more detailed analysis of the physics far below T c , Since we expect that So ~ Px~ l ^ ft we may 
expand the free energy © in j3J/Ti. Using Eqs. © we eliminate E and obtain U{B) = — tr(S 3 )/I27r/3 3 , resembling 
the SK-model where U(B) ~ —tr{B 2 )/j3 . The exponent reflects the spatial dimension D = 3 and is responsible for 
the shape of the pseudogap (p(E) ~ E ^ 1 ). In order to derive this result, it is more convenient to keep the self-energy 
in the formalism. Let us suppose that the replica symmetry is broken at the level of K steps. We represent the Parisi 
matrices as E = — Eo + J\_i Efei?,„ fc , where R mk consist of blocks of size mk on the diagonal with all entries equal 
to 1. Let us focus on the set C of the mi spins corresponding to one of the innermost blocks. These spins experience 
an effective field y created by all other spins. We describe its thermal fluctuations by a distribution P(y), which in 
the RS case was a simple Gaussian (see Eq. (JJjJ). In the case of continuous replica symmetry breaking, P(y) can in 
principle be obtained by integration of Parisi's differential equation using the methods of Refs. 20]. Here, we will 
only exploit the fact that the Coulomb glass is in a marginally stable state (the Hessian d 2 F/dB 2 has a vanishing 
eigenvalue in the replicon mode SB characterized by SB aa = and 5BR mh — for all /c), which imposes the constraint 

, d2/P(y) [2cosh(^/2)r = .gr^Eo)-^ 1 ^)- (8) 

Further, the innermost component of Eqs. © reads 

X = P ^ - (s a s^) a#£C =00i (So) 




FIG. 2: The distributions P(h) and P(y) of the local and thermodynamic fields, respectively. Since the latter allow for 
relaxation of the environment, the gap is much narrower. 
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(9) 



where we have introduced the charge susceptibility x, the analog of the zero field cooled susceptibility in spin glasses. 
Expanding g n for So//3 3> 1, one can see that at low temperatures these two equations only admit a solution if 
So ~ /3 and P(y) takes the scaling form 



T- 2 P{y = zT)^p{z) (T-0) 



(10) 



with p(z) ~ z 2 for j>1. This implies that the susceptibility obeys the scaling \ ^ T 2 , and the (static) screening 
length diverges at low temperatures as l sc — (47rx) -1 / 2 ~ T -1 . Note that % measures the charge response to a local 
potential change when the particles on other sites are allowed to readjust to the induced charge. Thus, it is associated 
with the thermodynamic local fields yt defined by (sj) = m, = tanh(/3yj/2)/2. While we expect \ t° control the 
hopping conductivity, tunneling experiments 5] probe the system on very short time scales, sampling the distribution 
P{h) of instantaneous local fields hi = ^ ■ JijSj. The thermal average of these local fields, (hi) = J^j <^ij m ji i s related 
to the thermodynamic field j/j via a Thouless-Anderson-Palmcr (TAP) equation 21], (ft,,) = y± + (s,-) /iq, where the 
Onsager term 
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(11) 



accounts for the extra polarizations induced by the presence of the charge (sj). For the consistency with the locator 
approximation, we have retained only terms corresponding to a local self-energy The deviation of the local field hi 
from its mean (hi) is essentially a Gaussian variable with width ho- More precisely, the relation 



P(h) = / dyP(y) 



cosh(/3h/2) e -0(h-y) 2 /2h o 



cos 



HPy/2) ^2Trho//3e-^ h o/^ 



(12) 



holds [22| , which generalizes a known result for the SK- model [23j . The tunneling density of states at zero bias then 
follows from Vq = [3 f dhP(h) [2 cosh(/3/i/2)]~ 2 . Eq. (|T2"|) implies that P(h) obeys a scaling analogous to Eq. JTDJ, 
and hence isq ~ T 2 . Generally, in order to make quantitative predictions, one needs to know the functional form of 
the field distributions. It turns out, however, that certain parameters arc not very sensitive to their details. It is 
convenient to assume a simple form P((h)) = a((h) + 7T 2 ) for the distribution of average fields, obtain P(y) via the 
TAP-equations and solve Eqs. I|8I9|) . This yields %, Vq and a as slowly varying functions of 7 [2J|: a ~ 0.204 — 0.00677, 
X w (22.27-0.8l7) T 2 , i/ « (2.178-0.0087)T 2 . The tunneling DOS i/ is roughly an order of magnitude smaller than 
the full susceptibility x, as is also evident from the typical distributions shown in Fig. [5] This agrees well with the 
experimental observation [12j that the susceptibilities governing tunneling and hopping transport differ significantly. 
The value of a should be compared to the Efros-Shklovskii prediction cues = 3/n sw 0.95 |3| which is larger than 
our estimate because their self-consistency argument imposes stability only with respect to single electron hops. By 
contrast, our estimate includes multiparticle constraints that decrease a below oles i n agreement with large-scale 
numerical simulations [25j. 

In conclusion, we have developed the locator approximation for Coulomb glasses, allowing us to include multiparticle 
correlations. We have used this formalism to provide evidence for a continuous glass transition below which the 



Coulomb glass gets stuck in a marginally stable state, resulting in subexponential relaxation dynamics and giving 
rise to the Coulomb gap. A priori, the locator approximation is justified for large disorder. However, as long as 
crystallization is prevented, a structural glass transition will provide sufficient self-generated disorder, so that we 
expect our results to hold at low temperatures even in the case of weak external disorder. We verified [22J that in this 
limit the local observables are still determined by the contribution from large scales and reveal the Efros-Shklovskii 
gap. Further, wc found that the locator approximation gives a significant decrease of the DOS with temperature 
already above T c , in agreement with numerics. Moreover, it predicts a discontinuous glass transition at a scale of 
T c ps 0.030 which depends, however, on the details of the cutoff at small scales. The validity of this prediction remains 
thus unclear. 

The locator approximation not only provides new insight into classical Coulomb glasses, but also allows for quan- 
titatively new predictions that go beyond the single particle theory. It thus sets the stage for further theoretical 
developments to understand the puzzles of correlated transport and glassy relaxation of these systems. For instance, 
it allows one to study the collective modes of the electrons which induce fluctuations in the local electric fields [llj 
and thus enhance the probability for resonant tunneling. Experiments indicate that such a mechanism might provide 
an alternative to phonon assisted tunneling, in particular at low temperatures when phonons freeze out [26j . Finally, 
extensions of the formalism to include quantum effects may be envisioned. 
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APPENDIX: JUSTIFICATION OF THE LOCATOR APPROXIMATION 



In order to justify the locator approximation we calculate the leading terms ignored by the locator approximation 
and show that their effect on the physical results is small in the limit W>1. To calculate the higher order terms we 
need a diagram technique. This is not convenient in the spin representation of the original problem 



H [{Si}] = 9 E SiJijSj + E ■ 



(13) 



but turns out to be rather simple in terms of the fields <; 
perform a Hubbard-Stratonovich transformation: 



conjugated to the spin variables Si that appear after we 



{Si} 



II' 
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|E, ij * i (/sj) ij Vj+E,(^+*) s -. 



(14) 



The field if3~ 1 cj) i can be thought of as the mean Coulomb potential at site i created by all other electrons |2(. 

In order to get rid of the disorder and to restore translational invariance, we apply the replica trick and calculate 
(Z n ) e , where () e denotes the average over the random energies e*. We assume the latter to be independent local 
variables with distribution -P(e). The replicated partition function, summed over spin degrees of freedom, reads 



(Z n ) t = / /'jJd0?e-^ w ..*fWJ)« 1 *?+E, 1 .McoBh(^4^ 



)] 



We expand the local terms ln[cosh( 2 )] with respect to <j>* and evaluate the disorder average by a cumulant 
expansion, 
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where -i—, C, ?7, A denote vertex coefficients: 
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[2cosh(/3e/2)] 2 
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6V27r/3iy 

Here we assumed a symmetric distribution of e, and dropped the zeroth order term in (|15|l since it is irrelevant in 
the replica limit n — > 0. The last equality in each equation gives the leading order in 1//3W for Gaussian disorder. 
Note that typically vertex coefficients scale as ~ 1/({3W), except for the coefficients c m corresponding to the replica 
conserving vertices (4>j) m - The latter are suppressed by additional powers of /3W, c m ~ l/(/3V7) m_1 because they 
correspond to disorder averages of exact derivatives of 1/cosh [/3e/2]. To establish the connection with the locator 
approximation, notice that in the leading order in 1//3W the quantity W coincides with the diagonal part of the self 
energy So of the replica symmetric theory. 

A systematic diagrammtic expansion is now obtained by taking [Gg -1 ]" 6 = < (PJ)^ 1 + Sij/fiW > 5 a + C,8ijl ab as 
bare propagator and treating all other terms in the cumulant expansion as interactions. In Fourier space, we have 



G a b (k) 



g (k)6 ao + 9l T 
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(16) 



9a(k) 

9i(k) = -CSo(fc) 



in 



The leading non-local contribution, T, 1 - 1 ' , to the replica diagonal part of the self-energy is a tripled propagator line, 
connecting two replica-conserving 4-vertices. Its fc-dependent part evaluates to 



E&(*)- s #(0) = 6 abV 2 j g§(r)(exp(ikr)-l)d 3 r 
= V 2 • 47r/? 3 <U \ 1 
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where r = (W /Air) 1 ' 2 . The local part is of the order of 
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As we will show below, the glass physics is dominated by long scales, k 
approximate the momentum dependent part of the self energy by 

-in 



< 



1/ro- For these momenta one can 



*$(*) - £i>) « -g^ 2 /? 3 (rofc) 2 + O ((r fc) 4 ) 



which results in an insignificant renormalization of the Coulomb interaction at temperatures of the order of the glass 
transition (/3 C ~ W 1 ^ 2 ). In this regime we can thus safely approximate the self-energy by its local part S^^(O). 

At very low temperatures the density of states at low energies is suppressed, this reduces the screening of the 
electric field and enhances the interactions. In this regime, the renormalizationsjjiight become more important. 
To estimate this effect we replace the screening part of the bare propagator, 1/W — /3([2cosh(/3e/2] _2 ) e by the 
average thermodynamic susceptibility \ — P J dyP(y)[2cosh(Py/2}~ 2 ~ 1//3 . This gives rise to the bare propagator 
<7o(fc) = P / (k 2 / in + x) with the longer screening length r$ — > (inx)^ 1 ^ 2 ~ /3. Further, as a consequence of the broken 
replica symmetry, the replica mixing vertices might also contribute to the renormalization of the replica diagonal part 
of the Greens function. At low temperatures the vertex coefficients generically scale as 1//3 . This might increase the 
importance of the non-local corrections and even make them marginally relevant below some temperature T*. By 
continuity we must have T* <T C . We do not know any physical argument that supports the appearance of the second 
temperature (energy) scale and it seems unlikclt to us that it happens. However, in order to check that non-local 
contributions remain parametrically small at all temperatures one needs to calculate them against the background of 
the full replica symmetry breaking solution. 

Incorporating the local term Sij6 a /W of the bare propagator into the self-energy S^ we get the full Greens function 



* ^3 



= \G~X = 



1 



-i ab 



[fi 'jy 1 - s 



(18) 



Before discussing the physical interpretation of this quantity in more detail we mention that the spin-spin correlation 
function follows easily via partial integration in l|14fl with respect to the conjugated fields, 
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(19) 



where £ = E" 1 . 

The locality of the self-energy suggests the mapping to an effective single site model, for which the same self-energy 
is assumed. For self-consistency we then have to require that the average local spin-spin correlators be the same. 
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(20) 



The equation l|19|) for spin-spin correlator is a central element in the mapping to a single site-model. Physically, 
this expression and the locator approximation in general, is based on the assumption that a typical spin interacts 
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FIG. 3: Leading diagrams for the computation of the four point function (</> i 0,-)'' 



with many thermally active neighbors. The calculations above show that this assumption is satisfied by a typical 
spin that contributes to the density of states and other thermodynamic properties. This, however, does not exclude 
the presence of ther types of low-energy excitations. For instance, an occupied and a nearby empty site may form a 
strongly coupled dipole, in which the electron can hop to the empty site at very low energy expense. Such pairs of sites 
are strongly coupled, and are not described by the locator approximation. The small effect of non-local corrections 
evaluated here physically mean that such dipoles renormalize weakly the dielectric susceptibility without interfering 
strongly in thermodynamics of carrier sites. The effect of these dipoles on other properties remain an open question. 



Fluctuating correlations and criticality 



The replica diagonal part of the Green function (|21(l is simply the disorder average (or equivalently, the average 
over sites at fixed distance) of the connected correlation function. Assuming a local self-energy, we explicitly obtain 



(^<t>i)=(Mj) c = 0- 



,—r/ro 



(21) 



where the overbar denotes the disorder average, and r = y / So/(4'?r/3) (= (W/iir) 1 ^ 2 in the limit of large W and at 
temperatures above T c ). This expression may be misleading since it seems to suggest that potentials are screened on 
average on a distance rg. However, we have to bear in mind that i|21|) describes only the disorder averaged correlations. 



To study the fluctuations around this average, we have to calculate 
the first term can be calculated as 



i'f'j) ~ i't'i't'j) ■ Within the replica formalism 
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(22) 



a,b,c and d being all distinct replicas that see the same disorder configuration. (In the glass phase one has to choose 
all indices in the same innermost block, which imposes that the two replicas belong to the same metastable state.) 
This four point function can be evaluated by the diagram technique, the leading diagrams being shown in Fig. 

The correlation function 122JI is given by the coefficient, C, of 8 ac 5bd of the full expression for 
easily be calculated by resumming a geometric series, which yields 






It can 



C(r) = / d 3 re 2fcr (0(O)0(r)) 2 = £ (2) (fc) 
Here, T,( 2 \k) denotes the insertion of a doubled propagator line, 



+ AS (2) (fc) + (AS (2) (/fc)) 2 + 



1/E( 2 )(fc)-A' 



(23) 



S (2) W = J 9 2 o(r)e lkr d 3 r = 2nr (l - M! + O((r fc) 4 ) 



Using the expressions at high temperature in the large VF-limit, one easily verifies that at the critical temperature 
T c = I / F _1 / 2 /[6(2/7r) 1 / 4 ] (cf. Eq. ©), the fc-independent terms in (|2*5)l cancel, leading to a power law decay of the 
potential correlations, 



<(/rV(o))(/rV(0)) r=T ~~-^ 



1 J3 C 
W 2 i 



(24) 



We separated a prefactor l/W 2 which has the interpretation of the probability that both sites at and r are thermally 
active. 

Note that the critical temperature is exactly the same as we found by mapping to a single-site problem. In the 
original lattice, the replicon instability receives a natural interpretation: The correlations in the Coulomb potential 
created by the electron configuration become critical at T c where they only decay with a power law. Below T c , the 
phase space splits into an exponential number of metastable states (ergodic components) each of which is characterized 
by the finite expectation values ipi of the conjugated potential fields ipi — iff {4>i)- From Eq. (|I4|I we see that 
within such a metastable state we should consider a + ^ as the effective field on the site i to be used to calculate 
the vertex coefficients. Anticipating an Efros-Shklovskii-type distribution for those effective fields then suggests that 
the vertex coefficients averaged over sites, scale as 1//3 as mentioned above. In particular, the mass term of the 
propagator will be replaced by the average susceptibility corresponding to the distribution of effective fields e^ + Vv 

At this point the diagram technique in the original lattice becomes very complicated since translational invariance 
is spontaneously broken by the emergence of spontaneous expectation values. The advantage of the mapping to a 
single site problem is now obvious: The replica structure implicitly takes the statistics of a spin's environment into 
account, and allows, e.g., for a self-consistent determination of the distribution of local fields. 

Spin-spin correlations and inhomogeneous charge response 

The criticality of the correlations of potential fields translates directly to the criticality of spins. The fact that the 
glass phase is marginally stable has the simple interpretation that these correlations remain long range (power law) 
throughout the low temperature phase. This self-organized criticality has an important consequence: One can show 
that by imposing the spin on site i to take a definite value one induces a polarization at site j which is proportional 
to (si-Sj) c . The above results imply that this response of the system in the glass phase decays only as a power law 
with distance, i.e., screening is almost absent. Furthermore, we note that the polarization induced on different sites 
do not have a definite sign. It should be clear from these considerations that typical spins interact with a large 
number of effective neighbors, which in turn self-consistently confirms the initial assumption suggesting the locator 
approximation. 

APPENDIX: THE REPLICON INSTABILITY 

Within the effective single-site model, the glass transition © and the marginal stability condition (JSJ both derive 
from the vanishing of the eigenvalue of the Hessian d 2 F/dB 2 in the replicon mode. The variation of the spin-part in 
(@J is standard and yields 



f°° 1 

S 2 p Fspm = -tr(SB 2 )j_JyP(y) [2cosHM2)]A . 



(25) 



When varying U(B), 10, we need to take into account that E is a function of B. The sclfconsistency condition 

f 



B-H 



= Si(-S) (26) 



imposes, upon variation in the replicon mode, that (SB — (5£)g 2 (Eo) — — <72(So)#£. The variation of U(B) thus 
evaluates to 

5 2 U(B) = -g 2 (Z )tr(SB - effi) 2 + . g2 (£ )ir(<ffi) 2 = ^^l, (27) 

01 ( s o) - .9 2 ( s o) 

which together with (|2*5|) leads to 10 and (JSJ) . 

APPENDIX: GENERALIZED ONSAGER TERM 

We have shown above that the spin so at site polarizes its environment in a large spatial region. As inis well- 
known from the TAP-approach to spin glasses, in order to obtain the thermodynamic field j/q = 2(3~ tanh~ (2mo), 



10 



the back reaction of this polarization on the spin itself has to be subtracted from the thermally averaged field 
(ho) =-T,jjto^ojmj, 



2/o = (ho) 



s h o , 



(28) 



where ho is the famous Onsager back reaction. The usual term ho = "YljJijXji familiar from spin glasses, has to 
be generalized to the case of Coulomb glasses where this expression is clearly divergent. Indeed, to obtain a finite 
response we have to sum up all higher order polarizations, 



hO — / j ^JijiXj-i <yjii / j •JijiXj l JjihXj 2 '-'k 



(29) 



n 



nj2 



their alternating sign reflecting the antiferromagnetic nature of the Coulomb interactions [2j. Approximating the 
on-site susceptibilities \j by their average \ (which is justified since we average over a very large number of spins), 
we may perform the sum 



ho = Tr- 



J 2 



x-' + J 



j3Tr 



J 2 



PJ + Z a 



27rVSn ~ T, 



(30) 



The last approximations are valid at low temperatures. 

The distribution of instantaneous fields can be obtained from 



V^J 2tt 



(31) 



We have only retained the first two cumulants, to be consistent within the locator approximation. From the generalized 
TAP-equations l|28|l we may identify the first cumulant as 



( h i) Si=s =Vi + sho, 
The second cumulant is almost insensitive to the value of the spin at site z, and evaluates to 



"YjJij (SjS k ) c Jki = 
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(32) 



(33) 



in the locator approximation. 

To carry out the site average, we have to weight the pairs (y, s), according to their joint probability density 
P( 2 /)exp(/3y S )/(2cosh(/3y/2), 



P(M) = JdyP(y) 

Performing the A-integral we find 
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dX 



e sf3y e i\(h-y-sH(!3)) e -^H(p)/2l3 



P(h, a) = / dyP(y) 



2cosh(f3y/2) J 2tt 

e s P y exp[-/3(/i -y- sh ) 2 /2h 



2cosh(/3y/2) 
and summing over s we find the local field distribution 
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(34) 



(35) 



,,.,. / . t>{ ,cosh(f3h/2)exp[-P(h-y) 2 /2h ] , n , lQ , ,„.. 

P{H) = I dvP{v) cosHfly/2) ^hoTW 2 6Xp(/3/lo/8) - (36) 



